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Abstract 

Spatio-temporally chaotic dynamics of a classical field can be described by means 
of an infinite hierarchy of its unstable spatio-temporally periodic solutions. The 
periodic orbit theory yields the global averages characterizing the chaotic dynamics, 
as well as the starting semiclassical approximation to the quantum theory. 

New methods for computing corrections to the semiclassical approximation are 
developed; in particular, a nonlinear field transformation yields the perturbative 
corrections in a form more compact than the Feynman diagram expansions. 
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Formulated in 1946-49 and tested through 1970's, quantum electrodynamics 
takes free electrons and photons as its point of departure, with nonlinear effects 
taken in account perturbatively in terms of Feynman diagrams, as corrections 
of order (a/vr)" = (0.002322819 . . .)". QED is a wildly successful theory, with 
Kinoshita's [1] calculation of the electron magnetic moment 




agreeing with Dehmelt's experiments [2] to 12 significant digits. 

Quantum chromodynamics perturbative calculations seemed the natural next 
step, the only new feature being the gluon-gluon interactions. However, in this 
case the Feynman-diagrammatic expansions for observables such as the meson 
and hadron masses 



(observable) = ^ («qcd)" 
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failed us utterly, perhaps because the expansion parameter is of order 1. I 
say perhaps, because more likely the error in this case is thinking in terms 
of quarks and gluons in the first place. Strongly nonlinear field theories re- 
quire radically different approaches, and in 1970's, with a deeper appreciation 
of the connections between field theory and statistical mechanics, their re- 
examination led to path integral formulations such as the lattice QCD [3]. In 
lattice theories quantum fiuctuations explore the full gauge group manifold, 
and classical dynamics of Yang-Mills fields plays no role. 

We propose to re-examine here the path integral formulation and the role that 
the classical solutions play in quantization of strongly nonlinear fields. In the 
path integral formulation of a field theory the dominant contributions come 
from saddlepoints, the classical solutions of equations of motion. Usually one 
imagines one dominant saddle point, the "vacuum": 



The Feynman diagrams of QED and QCD are nothing more than a scheme 
to compute the correction terms to this starting semiclassical, Gaussian sad- 
dlepoint approximation. But there might be other saddles. That field theories 
might have a rich repertoire of classical solutions became apparent with the 
discovery of instantons [4], analytic solutions of the classical SU{2) Yang- 
Mills equations of motion, and the realization that the associated instanton 
vacua receive contributions from countable oo's of saddles. What is not clear is 
whether these are the important classical saddles. Could it be that the strongly 
nonlinear theories are dominated by altogether different classical solutions? 

The search for the classical solutions of nonhnear field theories such as the 
Yang-Mills and gravity has so far been neither very successful nor very sys- 
tematic. In modern field theories the main emphasis has been on symmetries 
as guiding principles in writing down the actions. But writing down a differ- 
ential equation is only the start of the story; even for systems as simple as 
3 coupled ordinary differential equations one in general has no clue what the 
nature of the long time solutions might be. 

These are hard problems, and in explorations of modern field theories the dy- 
namics tends to be is neglected, and understandably so, because the wealth 
of the classical solutions of nonlinear systems can be truly bewildering. If 






one dominant extremum 



an infinity of instanton saddles 
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the classical behavior of these theories is anything like that of the field the- 
ories that describe the classical world — the hydrodynamics, the magneto- 
hydrodynamics, the Ginzburg-Landau system — there should be very many 
solutions, with very few of the important ones analytical in form; the strongly 
nonlinear classical field theories are turbulent, after all. Furthermore, there is 
not a dimmest hope that such solutions arc either beautiful or analytic, and 
there is not much enthusiasm for grinding out numerical solutions as long as 
one lacks ideas as what to do with them. 

By late 1970's it was generally understood that even the simplest nonlinear 

systems exhibit chaos. Chaos is the norm also for generic Hamiltonian flows, 
and for path integrals that implies that instead of a few, or countably few 
saddles, classical solutions populate fractal sets of saddles. 



For the path-integral formulation of quantum mechanics such solutions were 
discovered and accounted for by Gutzwiller [5] in late 1960's. In this framework 
the spectrum of the theory is computed from a set of its unstable classical 
periodic solutions. The new aspect is that the individual saddles for classically 
chaotic systems are nothing like the harmonic oscillator degrees of freedom, 
the quarks and gluons of QCD — they are all unstable and highly nontrivial, 
accessible only by numerical techniques. 

So, if one is to develop a semiclassical field theory of systems that are classically 
chaotic or "turbulent" , the problem one faces is twofold 

(1) determine, classify, and order by relative importance the classical solu- 
tions of nonlinear field theories. 

(2) develop methods for calculating perturbative corrections to the corre- 
sponding classical saddles. 

Our purpose here is to give an overview over the status of this program — for 
details the reader is referred to the literature cited. 

The first task, a systematic exploration of solutions of field theory has so 
far been implemented only for one of the very simplest field theories, the 
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1-dimensional Kuramoto-Sivashinsky system. Wc sketch below how its spatio- 
temporally chaotic dynamics can be described in terms of spatio-temporally 
recurrent unstable patterns. 

For the second task, the theory of perturbative corrections, we shall turn to 
an even simpler system; a weakly stochastic mapping in 1-dimension. The new 
aspect of the theory is that now the corrections have to be computed saddle 
by saddle. In sect. 3 to sect. 6 we discuss three distinct methods for their 
evaluation. 



1 Unstable recurrent patterns in classical field theories 

Field theories such as 4-dimensional QCD or gravity have many dimensions, 
symmetries, tensorial indices. They are far too complicated for exploratory 
forays into this forbidding terrain. We start instead by taking a simple spatio- 
temporally chaotic nonlinear system of physical interest, and investigate the 
nature of its solutions. 

One of the simplest and extensively studied spatially extended dynamical sys- 
tems is the Kuramoto-Sivashinsky system [6] 



which arises as an amplitude equation for interfacial instabilities in a variety 
of contexts. The "flame front" u{x,t) has compact support, with x e [0, 27r] a 
periodic space coordinate. The term makes this a nonlinear system, t is the 
time, and v is a fourth-order "viscosity" damping parameter that irons out 
any sharp features. Numerical simulations demonstrate that as the viscosity 
decreases (or the size of the system increases), the "flame front" becomes in- 
creasingly unstable and turbulent. The task of the theory is to describe this 
spatio-temporal turbulence and yield quantitative predictions for its measur- 
able consequences. 

Armed with a computer and a great deal of skill, one can obtain a numerical 
solution to a nonlinear PDE. The real question is; once a solution is found, 
what is to be done with it? The periodic orbit theory is an answer to this 
question. 

Dynamics drives a given spatially extended system through a repertoire of 
unstable patterns; as we watch a "turbulent" system evolve, every so often we 
catch a glimpse of a familiar pattern: 
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For any finite spatial resolution, the system follows approximately for a finite 
time a pattern belonging to a finite alphabet of admissible patterns, and the 
long term dynamics can be thought of as a walk through the space of such 
patterns, just as chaotic dynamics with a low dimensional attractor can be 
thought of as a succession of nearly periodic (but unstable) motions. The 
periodic orbit provides the machinery that converts this intuitive picture into 
precise calculation scheme that extracts asymptotic time predictions from the 
short time dynamics. For extended systems the theory gives a description of 
the asymptotics of partial differential equations in terms of recurrent spatio- 
temporal patterns. 

Putkaradze has proposed that the Kuramoto-Sivashinsky system (1) be used as 
a laboratory for exploring such ideas. We now summarize the results obtained 
so far in this direction by Christiansen et al. [7] and Zoldi and Greenside [8]. 

The solution u{x, t) — u{x + 27r, t) is periodic on the x e [0, 27r] interval, so 
one (but by no means only) way to solve such equations is to expand u{x, t) 
in a discrete spatial Fourier series 

+ 00 

u{x,t)=t ^ afe(t)e^^^ (2) 

fe=— oo 

Restrict the consideration to the subspace of odd solutions u{x, t) = —u{—x, t) 
for which Ofe are real. Substitution of (2) into (1) yields the infinite ladder of 
evolution equations for the Fourier coefficients a^: 

oo 

Ofc = {k'^ - i^k'^)ak -k Yl (^mak-m ■ (3) 

m=— oo 

u{x, i) = is a fixed point of (1), with the k'^u < 1 long wavelength modes of 
this fixed point linearly unstable, and the short wavelength modes stable. For 
1/ > 1, u{x,t) = is the globally attractive stable fixed point; starting with 
V = 1 the solutions go through a rich sequence of bifurcations, and myriad 
unstable periodic solutions whose number grows exponentially with time. 

The essential limitation on the numerical studies undertaken so far have been 
computational constraints: in truncation of high modes in the expansion (3), 
sufficiently many have to be retained to ensure the dynamics is accurately rep- 
resented. Christiansen et al. [7] have examined the dynamics for values of the 
damping parameter close to the onset of chaos, while Zoldi and Greenside [8] 
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Fig. 1. (a) The return map Sn+i = /(sn) constructed from periodic solutions of the 
Kuramoto-Sivashinsky equations (1), = 0.029910, with s the distance measured 
along the unstable manifold of the fixed point 1. Periodic points and 01 are also 
indicated, (b) Lyapunov exponents Xk versus k for the periodic orbit 1 compared 
with the stability eigenvalues of the u{x, t) = stationary solution k'^ — uk"^. for 
A: > 8 lie below the numerical accuracy of integration and are not meaningful. From 
ref. [7]. 

have explored somewhat more turbulent values of v. With improvement of nu- 
merical codes considerably more turbulent regimes should become accessible. 

One pleasant surprise is that even though one is dealing with (infinite dimen- 
sional) PDEs, for these strong dissipation values of parameters the spatio- 
temporal chaos is sufficiently weak that the flow can be visualised as an ap- 
proximately 1-dimensional Poincare return map s — > f{s) from the unstable 
manifold of the shortest periodic point onto its neighborhood, see figure 1(a). 
This representation makes it possible to systematically determine all nearby 
periodic solutions up to a given maximal period. 

So far some 1,000 prime cycles have been determined numerically for various 
values of viscosity. The rapid contraction in the nonleading eigendirections is 
illustrated in figure 1(b) by the plot of the first 16 eigenvalues of the 1-cycle. 
As the length of the orbit increases, the magnitude of contracting eigenvalues 
falls off very quickly. In figure 2 we plot Uq{x, t) corresponding to the 0-cycle. 
The difference between this solution and the other shortest period solution is 
of the order of 50% of a typical variation in the amplitude of u{x,t), so the 
chaotic dynamics is already exploring a sizable swath in the space of possible 
patterns even so close to the onset of spatio-temporal chaos. Other solutions, 
plotted in the configuration space, exhibit the same overall gross structure. 
Together they form the repertoire of the recurrent spatio-temporal patterns 
that is being explored by the turbulent dynamics. 
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(a) 




Fig. 2. (a) Spatio-temporally periodic solution uo{x, t) of the Kuramoto-Sivashinsky 
system, viscosity parameter u = 0.029910. (b) The difference between the two 
shortest period spatio-temporally periodic solutions uo{x,tTo) and ui{x,tTi). Prom 



2 Periodic orbit theory 

Now we turn to the central issue; qualitatively, these solutions demonstrate 
that the recurrent patterns program can be implemented, but how is this in- 
formation to be used quantitatively? This is what the periodic orbit theory is 
about; it offers the machinery that assembles the topological and the quanti- 
tative information about individual solutions into accurate predictions about 
measurable global averages, such as the Lyapunov exponents and correlation 
functions. 

Very briefly (for a detailed exposition the reader is referred to ref. [9]), the 
task of any theory that aspires to be a theory of chaotic, turbulent systems is 
is to predict the value of an "observable" a from the spatial and time averages 
evaluated along dynamical trajectories x{t) 



The key idea of the periodic orbit theory is to extract this average from the 
leading eigenvalue of the evolution operator 



ref. [7]. 








via the trace formula 



tr 
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EE 



oo 



T,S{t-rTp) 



,r/3A, 



(4) 



p 



p r=l 



det (i - J;) 



p 
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which relates the spectrum of the evolution operator to a sum over prime 
periodic solutions p of the dynamical system and their repeats r. 



What does this formula mean? Prime cycles partition the dynamical space into 
neighborhoods, each cycle enclosed by a tube whose volume is the product of 
its length Tp and its thickness |dct(l — Jp)|~^. The trace picks up a periodic 
orbit contribution only when the time t equals a prime period or its repeat, a 
constraint enforced here by 6{t — rTp). 3p is the linear stability of cycle p, so 



for long cycles 
contribution of 



det (l — ~ (product of expanding eigenvalues), and the 
ong and very unstable cycles are exponentially small compared 
to the short cycles which dominate trace formulas. The number of contracting 
directions and the overall dimension of the dynamical space is immaterial; 
that is why the theory can also be applied to PDEs. All this information is 
purely geometric, intrinsic to the flow, coordinate reparametrization invariant, 
and the same for any average one might wish to compute. The information 
related to a specific observable is carried by the weight e^'^p, the periodic orbit 
estimate of the contribution of e^^ from the p-cycle neighborhood. 

The intuitive meaning of a trace formula is that it expresses the average (e^"^^ 
as a discretized integral 




over the dynamical space partitioned topologically into a repertoire of spatio- 
temporal patterns, each weighted by the likelihood of pattern's occurrence in 
the long time evolution of the system. 



Periodic solutions are important because they form the skeleton of the invari- 
ant set of the long time dynamics, with cycles ordered hierarchically; short cy- 
cles give good approximations to the invariant set, longer cycles refinements. 
Errors due to neglecting long cycles can be bounded, and for nice hyperbolic 
systems they fall off exponentially or even super-exponentially with the cutoff 
cycle length [10]. Short cycles can be accurately determined and global aver- 
ages (such as Lyapunov exponents and escape rates) can be computed from 
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short cycles by means of cycle expansions. 

The Kuramoto-Sivashinsky periodic orbit calculations of Lyapunov exponents 
and escape rates [7] demonstrate that the periodic orbit theory predicts observ- 
able averages for deterministic but classically chaotic spatio-temporal systems. 
The main problem today is not how to compute such averages — periodic orbit 
theory as well as direct numerical simulations can handle that — but rather 
that there is no consensus on what the sensible experimental observables worth 
are predicting. 

It should be obvious, and it still needs to be said: the spatio-temporally pe- 
riodic solutions are not to be thought of as eigenmodes, a good linear basis 
for expressing solutions of the equations of motion. Something like a dilute in- 
stant approximation makes no sense at all for strongly nonlinear systems that 
we are considering here. As the equations are nonlinear, the periodic solutions 
are in no sense additive, and their linear superpositions are not solutions. 




7^ u(x,t) 



a solution a solution not a solution 

Instead, it is the trace formulas and spectral determinants of the periodic 
orbit theory that prescribe how the repertoire of admissible spatio-temporal 
patterns is to be systematically explored, and how these solutions are to be 
put together in order to predict measurable observables. 

Suppose that the above program is successfully carried out for classical so- 
lutions of some field theory. What are we to make of this information if we 
are interested in the quantum behavior of the system? In the semiclassical 
quantization the classical solutions are the starting approximation. 



3 Stochastic evolution 

For the same pragmatic reasons that we found it profitable to shy away from 
facing the 4-dimensional QCD head on in the above exploratory foray into 
a strongly nonlinear field theory, we shall start out by trying to understand 
the structure of perturbative corrections for systems radically simpler than a 
full-fiedgcd quantum field theory. First, instead of perturbative corrections to 
the quantum problem, we shall start by exploring the perturbative corrections 
to weakly stochastic fiows. Second, instead of continuous time fiows, we shall 
start by a study of a discrete time process. 
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For discrete time dynamics a Langevin trajectory in presence of additive noise 
is generated by iteration 



Xn+l = f{Xn) + CrCn , (5) 

where f{x) is a map, ^„ a random variable, and a parametrizes the noise 
strength. In what follows we assume that ^„ are uncorrelated, and that the 
mapping f(x) is one-dimensional and expanding, but we expect that the form 
of the results will remain the same for higher dimensions, including the field 
theory example of the preceding section. 

Tracking an individual noisy trajectory does not make much sense; what makes 
sense is the Fokker-Planck formulation, where one considers instead evolution 
of an ensemble of trajectories. An initial density of trajectories (l)o{x) evolves 
with time as 

0n+i(y) = (>C o 0„) (y) = j dxC{y, (6) 



where C is the evolution operator 

x)^ j 5{y - f(x) - a^Pim = (^''P - /(^))] , (7) 



and ^ random variable with the normalized distribution P{C), centered on 

If the noise is weak, the goal of the theory is to compute the perturbative 
corrections to the eigenvalues u of C order by order in the noise strength a, 



m=0 



One way to get at the spectrum of C is to consider the discrete Laplace 
transform of C^, or the resolvent 

V ;2"tr = tr = V (8) 



which has a pole at every z — v^^. 

The effects of weak noise are of interest in their own right, as any deterministic 
evolution that occurs in nature is affected by noise. However, what is most 
important in the present context is the fact that the form of perturbative 
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corrections for the stochastic problem is the same as for the quantum problem, 
and still the actual calculations are sufficiently simple that one can explore 
many more orders in perturbation theory than would be possible for a full- 
fledged field theory, and develop new perturbative methods. 

The first method we try is the standard Feynman-diagrammatic expansion. 
For semiclassical quantum mechanics of a classically chaotic system such cal- 
culation was first carried out by Gaspard [11]. The stochastic version described 
here, implemented by Dettmann [12], reveals features not so readily apparent 
in the quantum calculation. 

The Feynman diagram method becomes unwieldy at higher orders. The second 
method, introduced by Vattay [13], is based on Rugh's [10] explicit matrix 
representation of the evolution operator. If one is interested in evaluating 
numerically many orders of perturbation theory and many eigenvalues, this 
method is unsurpassed. 

The third approach, the smooth conjugations introduced by Mainieri [14], is 
perhaps an altogether new idea in field theory. In this approach the neighbor- 
hood of each saddlepoint is rectified by an appropriate nonlinear field transfor- 
mation, with the focus shifted from the dynamics in the original field variables 
to the properties of the conjugacy transformation. The expressions obtained 
are equivalent to sums of Feynman diagrams, but are more compact. 



4 Feynman diagrammatic expansions 



We start our computation of the weak noise corrections to the spectrum of L 
by calculating the trace of the n-th iterate of the stochastic evolution oper- 
ator L. A convenient choice of noise is Gaussian, P(^) = /^/-s/Stt , with the 
trace given by an n-dimensional integral on n points along a discrete periodic 
chain 



tr>C" = ^ dxo- ■ -dxri-iCixo^Xn-i) • • • C{xi,xo) 

1 

'2^ 



'A exp I [Xa+l - f{Xa)Y 



Xn^XQ, [dx] = n -^f^^ ■ (9) 

The choice of Gaussian noise is not essential, as the methods that we develop 
here apply equally well to other noise distributions, and more generally to 
the space dependent noise distributions P{x,^). As the neighborhood of any 
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trajectory is nonlincarly distorted by the flow, the integrated noise is anyway 
never Gaussian, but colored. 

If the classical dynamics is hyperbolic, periodic solutions of given flnite period 
n are isolated. Furthermore, if the noise broadening a is sufficiently weak they 
remain distinct, and the dominant contributions come from neighborhoods of 
periodic points, the tubes sketched in the trace formula (4). In the saddlepoint 
approximation the trace (9) is given by the sum over neighborhoods of periodic 
points 

oo 

t.C-^t.C\^^ Y: e^= = ^n,^5„,„^,e^^'-. (10) 
xcsFix/" P r=i 

As traces are cyclic, e^" is the same for all periodic points in a given cycle, 
independent of the choice of the starting point Xc, and the periodic point 
sum can be rewritten in terms of prime cycles p and their repeats. In the 
deterministic, cr — > limit this is the discrete time version of the classical 
trace formula (4) . Effects such as noise induced tunnelling are not included in 
the weak noise approximation. 

We now turn to the evaluation of Wpr, the weight of the r-th repeat of prime 
cycle p. The contribution of the cycle point Xa neighborhood is best expressed 
in an intrinsic coordinate system, by centering the coordinate system on the 
cycle points, 

Xa^ Xa + (f)a- (H) 

Prom now on Xa will refer to the position of the a-th periodic point, 0^ to the 

deviation of the noisy trajectory from the deterministic one, fa{4'a) to the map 
(5) centered on the a-th cycle point, and f^"^^ to its m-th derivative evaluated 
at the a-th cycle point: 

faM = f(Xa + K) - Xa+1 , fa = /'(^a), = /" (^a), " " " ■ (12) 

Rewriting the trace in vector notation, with x and f{x) n-dimensional column 

vectors with components Xa and /(x^) respectively, expanding / in Taylor 
series around each of the periodic points in the orbit of Xc, separating out the 
quadratic part and integrating we obtain 

e'^^ = J [d<P] e-(^- V-v^'W)'/2.^ ^ J ^"z^<^ + (-) 

c c 

= IdetAI J m eE^^(AV"W)'=e-^V2.^ ^13) 
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The [nxn] matrix A arises from the quadratic part of the exponent, while all 
higher powers of (pa are collected in V{(f)): 

oo rb"'-^^ 

Aa.V. = -/>a + , V{<P) = 7Z^, ' (14) 

a m=2 [m -\- I). 

The saddlepoint expansion is most conveniently evaluated in terms of Feynman 
diagrams, by drawing A as a directed line Aab = — *— , and the derivatives 
of V as the "interaction" vertices Q 

fl = fa = 

In the jargon of field theory, A is the "free propagator". Its determinant 

-1 n— 1 

|detA| = — — -, Ac = Ufa (15) 

M^c -'-I a=0 

is the 1-dimensional version of the classical stability weight |det(l — J)P^ in 
(4), with Ac the stability of the n-cycle going through the periodic point Xc- 

Standard methods [15] now yield the perturbation expansion in terms of the 
connected "vacuum bubbles" 

oo 

Wc = - In |Ac - 1| + E Wc,2k(r''' (16) 
fc=i 

In the usual field-theoretic calculations the Wc,o term corresponds to an over- 
all volume term that cancels out in the expectation values. In contrast, as 
explained in sect. 2, here the e^'='° = |Ac — 1|~^ term is the classical vol- 
ume of cycle c. Not only does this weight not cancel out in the expectation 
value formulas, it plays the key role both in classical and semiclassical trace 
formulas. 

In the diagrams sketched above a propagator line connects Xa at time a with 
Xh at later time 6 by a deterministic trajectory. At time b noise induces a 
kick whose strength depends on the local curvature of the flow. A penalty of a 

^ PC: fix birdtrack height 




13 



factor a is paid, m — 1 deterministic trajectories originate in the neighborhood 
of ,T5 from vertex ^'^'"^(xb), and the process repeats itself, each vertex carrying 
a penalty of cr, and higher derivatives of the f^. Summing over all noise kick 
sequences encoded by a given diagram and using the periodicity of the trace 
integral (9) Dettmann [12] obtains expressions such as 

=6+1 





This particular sum is the (^^^^ Feynman diagram correction to r-th 

repeat of prime cycle p. More algebra leads to similar contributions from the 
remaining diagrams. But the overall result is surprising; the dependence on 
the repeat number r factorizes, with each diagram yielding the same prefac- 
tor depending only on A^. This remarkable fact will be explained in sect. 6. 
The result of the Feynman-diagrammatic calculations is the stochastic trace 
formula 



tr 



1-zC 



_ nptp,k _ Z^^ ^^W+0(^4) 



where tp^k is the /c-th local eigenvalue evaluated on the p cycle. The determin- 
istic, o" = part of this formula is the stochastic equivalent of the Gutzwiller 
semiclassical trace formula [5]. The o"^ correction w^^\ is the stochastic ana- 
logue of Gaspard's h correction [11]. At the moment the explicit formula is 
sufficiently unenlightening that we postpone writing it down to sect. 6. 

While the diagrams are standard, the chaotic field theory calculations are 
considerably more demanding than is usually the case in field theory. Here 
there is no translational invariance along the chain, so the vertex strength 
depends on the position, and the free propagator is not diagonalized by a 
Fourier transform. Furthermore, here one is neither "quantizing" around a 
trivial vacuum, nor a countable infinity of analytically explicit soliton saddles, 
but around an infinity of nontrivial unstable hyperbolic saddles. 

Two aspects of the above perturbative results are a priori far from obvious: 
(a) that the structure of the periodic orbit theory should survive introduction 
of noise, and (b) a more subtle and surprising result, repeats of prime cycles 
can be re-summed and theory reduced to the dynamical zeta functions and 
spectral determinants of the same form as for deterministic systems. 

Pushing the Feynman-diagrammatic approach to higher orders is laborious, 
and has not been attempted for this class of problems. As we shall now see, it 
is not smart to keep pushing it, either, as one can compute many more orders 
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of perturbation theory by means of a matrix representation for L. 



5 Evolution operator in a matrix representation 



An expanding map f{x) takes an initial smooth distribution (\){x) defined on 
a subinterval, stretches it out and overlays it over a larger interval. Repetition 
of this process smoothes the initial distribution so it is natural to con- 

centrate on smooth distributions 0„(x), and represent them by their Taylor 
scries. By expanding both and (f)n+i{y) in (6) in Taylor series Rugh [10] 

derived a matrix representation of the evolution operator 

dxC{y,x)— = J2—iT'^m'm, m, m' = 0, 1, 2, . . . 
ml , m'\ 



which maps the x^ component of the density of trajectories in (6) to 

the y"* component of the density 0„+i(y) one time step later. The matrix 
elements follow by differentiating both sides with d^' jdy^' and evaluating 
the integral 

^m'm-g^JdxCiy,x)- (19) 



In (7) we have written the evolution operator C in terms of the Dirac delta 
function in order to emphasize that in the weak noise limit the stochastic 
trajectories are concentrated along the classical trajectory y — f{x). Hence it 
is natural to expand the kernel in a Taylor series [17] in a 

jC{y, x) = 5(y - f(x)) + ^ ^--^5^^\y - f{x)) / ^^(0^^ , (20) 

n=2 ^- 



where 5*'"^(y) = ^^S{y) . This yields a representation of the evolution operator 
centered along the classical trajectory, dominated by the deterministic Pcrron- 
Frobenius operator 6{y — f{x)), with corrections given by derivatives of delta 
functions weighted by moments of the noise distribution Pn — J -P(0^"^^- We 
again center the coordinate system on the cycle points as in (11), and also 
introduce a notation for the operator (7) centered on the Xa — > Xa+i segment 
of the classical trajectory 

>Ca(0a+l, 4>a) = jO.{Xa+l + 4>a+l, Xa + 4>a) ■ 
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The weak noise expansion (20) for the a-th segment operator is given by 



jCai<f>a+l, (t>a) = K<t>a+1 - fa{<t>a)) + E ^-^PJ^''\(t>a+l - fa{<t>a)) ■ (21) 



n=2 



As the evolution operator has a simple 5-function form, the local matrix rep- 
resentation of Ca centered on the Xa — > Xa+i segment of the deterministic 
trajectory can be evaluated recursively in terms of derivatives of the map /: 



(I^a)m'm = E ^ ^ (Ba)m'+n,m , n = max(m - m', 0) 

d05('"')(0„+i-/a(0))^4 (22) 

m! 



1 /rf 1 



Lm 



\fL\ \d<t> f'a{<t>) ) rn\ 

The matrix elements vanish for m' < m, so B is a lower triangular matrix. The 
diagonal and the successive off-diagonal matrix elements are easily evaluated 
iteratively by computer algebra 

1 (^+2)!/: 



mfT' V 2m!|/^|(/^) 



m+2 



For chaotic systems the map is expanding, |/^| > 1. Hence the diagonal terms 
drop off exponentially, as l/l/aP"*"^, the terms below the diagonal fall off even 
faster, and truncating to a finite matrix introduces only exponentially small 
errors. 



The trace formula (8) takes now a matrix form 



tr 



zC 



Vn.tr , (23) 



where Lp = L„pL2 • • • Li is the contribution of the p cycle. The subscript sc 
is a reminder that this is a saddlepoint or semiclassical approximation, valid 
as an asymptotic series in the limit of weak noise. Vattay [16] interprets the 
local matrix representation of the evolution operator as follows. The matrix 
identity log dct = tr log together with the trace formula (23) yields 



det(l-^£)L^ = ndet(l-^"''L,), (24) 
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so in the saddlepoint approximation the spectrum of the global evolution oper- 
ator C is pieced together from the local spectra computed cycle-by-cycle on 
neighborhoods of individual prime cycles with periodic boundary conditions. 
The meaning of the A;-th term in the trace formula (18) is now clear; it is 
the k-th eigenvalue of the local evolution operator restricted to the p-th cycle 
neighborhood. 

Using this matrix representation Palla and S0ndergaard [13] were able to com- 
pute corrections to order cr^^, a feat simply impossible along the Feynman- 
diagrammatic line of attack. In retrospect, the matrix representation method 
for solving the stochastic evolution is eminently sensible — after all, that is 
the way one solves a close relative to stochastic PDEs, the Schrodinger equa- 
tion. What is new is that the problem is being solved locally, periodic orbit 
by periodic orbit, by translation to coordinates intrinsic to the periodic orbit. 
It is this natural local basis that makes the matrix representation so simple. 

Mainieri [14] takes this observation one step further; as the dynamics is non- 
linear, why not search for a nonlinear coordinate transformation that makes 
the intrinsic coordinates as simple as possible? 



6 Smooth conjugacies 

This step injects into field theory a method standard in the construction of 
normal forms for bifurcations [18]. The idea is to perform a smooth nonlinear 

coordinate transformation x = h{y), f{x) = h{g{h~^{x))) that flattens out the 
vicinity of a fixed point and makes the map linear in an open neighborhood. 



The key idea of fiattening the neighborhood of a saddlepoint can be traced 
back to Poincare's celestial mechanics, and is perhaps not something that a 
field theorist would instinctively hark to as a method of computing pertur- 
bative corrections. This local rectification of a map can be implemented only 
for isolated non- degenerate fixed points (otherwise higher terms are required 
by the normal form expansion around the point), and only in finite neighbor- 
hoods, as the conjugating functions in general have finite radia of convergence. 

We proceed in two steps. First, substitution of the weak noise perturbative 
expansion of the evolution operator (21) into the trace centered on cycle c 



f{x) g{y) =J-y. 




an arbitrary coordinatization 



intrinsic, fiat coordinates 
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generates products of derivatives of ^-functions: 

irc\ = --- + j m {■ ■ ■ - M')) ^^"'K'P' - fa-M) •••} + • 



The integrals are evaluated as in (22), yielding recursive derivative formulas 
such as 



d l_ 

dx y'{x) 



y^f{x)-x. (25) 



y=o 



or n-point integrals, with derivatives distributed over n different 5-functions. 

Next we linearize the neighborhood of the a-th cycle point. For a 1- dimensional 
map f{x) with a fixed point /(O) = of stability A — f (0) , \A\ ^ 1 we search 
for a smooth conjugation h{x) such that: 

f{x) = h{Ah-\x)), h{0) = 0, /i'(0) = l. (26) 



In higher dimensions A is replaced by the Jacobian matrix J. For a periodic 
orbit each point around the cycle has a differently distorted neighborhood, 
with differing second and higher derivatives, so the conjugation function ha 
has to be computed point by point. 

An explicit expression for ha in terms of / is obtained by iterating around the 
whole cycle, and using the chain rule (15) for the cycle stability Ap 

f:^{<P)^ha{A,h-\<f>)), (27) 

so each ha is given by some combination of fa derivatives along the cycle. 
Expand f{x) and h{x) 

fix) ^Ax + x^f2 + xVa + ■ ■ ■ , h{y) + y^h2 + y^h^ + ... , 

and equate recursively coefficients in the functional equation h{Ay) — f{h{y)) 
expansion 

CXD 

h{Au) - Ah{u) = J2 fm {h{u)r . (28) 

n=2 
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This yields the expansion for the conjugation function h in terms of the map- 
ping / 

A(A-l)' A2(A-1)(A2-1)' ■ ^29) 



The periodic orbit conjugating functions ha are obtained in the same way from 
(27), with proviso that the cycle stability is not marginal, |Ap| 7^ 1. 

What is gained by replacing the perturbation expansion in terms of by 
still messier perturbation expansion for the conjugacy function hi Once the 
neighborhood of a fixed point is linearized, the conjugation formula for the 
repeats of the map 

r{x) = h{A^h-\x)) 



can be used to compute derivatives of a function composed with itself r times. 
The expansion for arbitrary number of repeats depends on the conjugacy 
function h{x) computed for a single repeat, and all the dependence on the 
repeat number is carried by polynomials in A*", a result that emerged as a 
surprise in the Feynman diagrammatic approach of sect. 4. The integrals such 
as (25) evaluated on the r-th repeat of prime cycle p 

y{x) = r'{x) - X (30) 



have a simple dependence on the conjugating function h 



1 1 

1 1 
4! dy^ y'(0) 



A'' (1 + A'' 
(A'- - 1)^ 

(A^-l)4 2 



5A 



2r 



8A'' 



(A'- - 1)4 



-/12/13 



(31) 



,^A2'- + A'- + 1 
(A'- - 1)4 



The evaluation of n-point integrals is more subtle [14] . The final result of all 
these calculations is that expressions of form (31) depend on the conjugation 
function determined from the iterated map, with the saddlepoint approxima- 
tion to the spectral determinant given by 

00 

det(i-^£,)L=nn(i-Vfc) 

P k=0 



in terms of local p-cycle eigenvalues 
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|Ap|Ap 

a 

accurate up to order a*. w^^\ w^^'' are also computed by Dettmann, but we de- 
sist from citing them here; the reader is referred to ref. [14]. What is remarkable 
about these results is their simplicity when expressed in terms of the conju- 
gation function h, as opposed to the Feynman diagram sums, in which each 
diagram contributes a sum like the one in (17), or worse. Furthermore, both 
the conjugation and the matrix approaches are easily automatized, as they 
require only recursive evaluation of derivatives, as opposed to the handcrafted 
Feynman diagrammar. 

Simple minded as they might seem, discrete stochastic processes arc a great 
laboratory for testing ideas that would otherwise be hard to test. Dettmann, 
Palla and S0ndergaard have used a 1-dimensional repeller of bounded non- 
linearity and complete binary symbolic dynamics to check numerically the 
above results, and computed the leading eigenvalue of JC by no less than five 
different methods. As anticipated by Rugh [10], the evolution operator eigen- 
values converge super-exponentially with the cycle length; addition of cycles 
of period (n+1) to the set of all cycles up to length n doubles the number 
of significant digits in the perturbative prediction. However, as the series is 
asymptotic, for realistic values of the noise strength summations beyond all 
orders are needed [19]. 



7 Summciry 



The periodic orbit theory approach to turbulence is to visualize turbulence 
as a sequence of near recurrences in a repertoire of unstable spatio-temporal 
patterns. The investigations of the Kuramoto-Sivashinsky system discussed 
above are first steps in the direction of implementing this program. So far, 
existence of a hierarchy of spatio-temp orally periodic solutions of spatially 
extended nonlinear system has been demonstrated, and the periodic orbit 
theory has been tested in evaluation of global averages for such system. The 
parameter ranges tested so far probe the weakest nontrivial "turbulence", and 
it is an open question to what extent the approach remains implementable as 
the system goes more turbulent. 

The most important lesson of this investigation is that the unstable spatio- 
temporally periodic solutions do explore systematically the repertoire of ad- 
missible spatio-temporal patterns, with the trace and spectral determinants 
formulas and their cycle expansions being the proper tools for extraction of 
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quantitative predictions from the periodic orbits data. 

We formulat next a semiclassical perturbation theory for stochastic trace for- 
mulas with support on infinitely many chaotic saddles. The central object of 
the periodic orbit theory, the trace of the evolution operator, is a discrete 
path integral, similar to those found in field theory and statistical mechan- 
ics. The weak noise perturbation theory, likewise, resembles perturbative field 
theory, and can be cast into the standard field-theoretic language of Feyn- 
man diagrams. However, we found out that both the matrix and the nonlinear 
conjugacy perturbative methods are superior to the standard approach. In 
contrast to previous perturbative expansions around vacua and instanton so- 
lutions, the location and local properties of each saddlepoint must be found 
numerically. 

The key idea in the new formulation of perturbation theory is this: Instead of 
separating the action into quadratic and "interaction" parts, one first performs 
a nonlinear field transformation which turns the saddle point into an exact 
quadratic form. The price one pays for this is the Jacobian of the nonlinear 
field transformation — but it turns out that the perturbation expansion of this 
Jacobian in terms of the conjugating function is order-by-order more compact 
than the Feynman-diagrammatic expansion. 
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